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POPULATION INVERSION CALCULATIONS USING NEAR RESONANT 

CHARGE EXCHANGE AS A PUMPING MECHANISM 

by Donald L. Chubb and James R. Rose 
Lewis Research Center 


ABSTRACT 

Near resonance charge exchange between ions of a large ionization 
potential gas such as helium or neon and vapors of metals such as zinc, 
cadmium, selenium, or tellurium has produced laser action in the metal 
ion gas. The present investigation is a theoretical study of the possi- 
bility of obtaining population inversions in near resonant charge exchange 
systems (Xe-Ca, Xe-Mg, Xe-Sr, Xe-Ba, Ar-Mg, N-Ca) . The analysis is an 
initial value problem that utilizes rate equations for the densities of 
relevant levels of the laser gas (Ca, Ba, Mg, or Sr) and an electron 
energy equation. Electron excitation rates are calculated using the Bohr- 
Thomson approximation for the cross section. Approximations to experi- 
mental values of the electron ionization cross section and the ion-atom 
charge exchange cross section are used. Preliminary results have been 
obtained for the Ca-Xe system and show that it is possible to obtain gains 
greater than 1C>14 m - 3 with inversion times up to 8xl0~^ second. A pos- 
sible charge exchange laser system using a MPD arc plasma accelerator is 
also described. 


I. INTRODUCTION 

McGowan and Stebbings (ref. 1) proposed that near-resonance charge 
exchange is a possible way of obtaining an inversion in a gas laser. Ex- 
perimentally it was found that charge exchange was a pumping mechanism 
for producing inversions in He-Zn (refs. 2, 3, and 4), He-Cd (ref. 5), 
He-Se (refs. 6 and 7), He-Te (ref. 8), and Ne-Te (ref. 8). Theoretically 
there are many systems where near-resonance charge exchange can be used 
to produce an inversion. 

The present theoretical study investigates the possibility of obtain- 
ing inversions in a plasma of metals such as Mg, Ca, Sr, or Ba by charge 
exchange with inert gas ions such as He, Ne, Xe, Ar, or N. The analysis 
uses rate equations to calculate the density of the various ion levels 
considered and an energy equation to calculate the electron temperature. 
Both charge exchange and electron collisional excitation of the various 
ion levels are considered. 

In the next section a charge exchange laser system is described. 
Following that the analysis is presented. Then the important collision 
cross sections are discussed. Finally results of the analysis for the 
Ca-Xe system are discussed and conclusions are drawn. 
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II. CHARGE EXCHANGE ION LASER 

Near resonance charge exchange is described by the following reac- 
tion. 


A + + B 


AO 


A + 


B + + 
x 


AE 


(la) 


The excited ion, B-^, corresponds to the upper laser level of the metal 
being considered. The ground state ion, A , corresponds to the inert gas 
ion. The quantity is the excitation rate for the B"^ state by charge 

exchange with A + . The quantity AE is given by 

E = I.P. A - (l.P. B + E.P. B +) (lb) 


where I.P. is ionization potential E.P. is excitation potential. If 
AE is small, then the reaction rate, JJ^q, f° r this process will be com- 
parable, to the resonance process where AE = 0. Such charge exchange 
rates are large compared to electron and radiative excitation rates and 
therefore make it possible to preferentially excite the (B + )^ state. 

There are mary possible choices for A and B that satisfy the require- 
ment AE -* 0. It is a reaction of the type described by equation (la) 
that produced inversions in the experiments of references 2 to 8. 

Table I lists some interesting possible systems. The first column 
lists the charge exchange reaction that populates the upper laser level 
of the metal ion. Column two lists the possible laser transition and col- 
umn three gives the laser wavelength (ref. 18). The fourth column. gives 
the value of the energy defect, AE, and the last column is the upper bound 
on the efficiency of the system. This efficiency is defined as the ratio 
of the energy of the laser photon, hc/X, to the ionization potential of 
the inert gas, I.P. (species A in eq. (la)). 

In addition to the charge exchange reaction of equation (la) , the 
following charge exchange reaction may also produce inversions. 

+ * AB+ + 

A + B -*■ . A + B ± + AE (2a) 


Here, 

E = I.P.g - (I.P. A + E.P. ) (2b) 

Again there are many systems where AE -*■ 0 so that a large reaction rate 
can be expected. 
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Referring to equation (lb) we see that the A species must have a 
large ionization potential (for example, an inert ga^) in order that 
AE ■+ 0 and equation (la) be a near resonant reaction. For. equation (2a) 
to be a hear resonant reaction, however, the A species must have a small 
ionization potential (for example, a metal). 

Previous laser action based on equation (la) has been obtained in 
discharge tubes running on He or Ne seeded with a metal (refs. 2 to 8). 
Inherently such devices establish plasma conditions (electron temperature 
and density) necessary to maintain the discharge. Usually these same con- 
ditions are not proper for obtaining large gains. Hopefully, as power is 
increased to the discharge the gain should increase. However, as power is 
increased the gain increases to a maximum and then decreases (ref; 7). The 
reason for the decrease in gain with increasing power is that the electron 
temperature increases with power and finally becomes large enough so that 
electron collisions drive all the excited states toward their equilibrium 
values thus eliminating the inversion. 

The ideal system for a charge exchange laser would be one where .the 
plasma conditions for obtaining gain are independent of those necessary 
for maintaining the source of ions. One device that has the possibility 
of providing this independent control is a plasma accelerator called an MPD 
arc (refs. 9 and 10). The MPD arc produces a stream of high velocity 
(1.5 X 10^ m/sec) ions at densities up to 10^0 m - 3 ( re f. 11). A proposed 
system for obtaining gain using charge exchange with an MPD arc is illus- 
trated in figure 1. A seed metal is added to the exhaust of the arc. An 
inversion in the metal ion gas should be produced by means of the charge 
exchange reaction described by equation (la) . Since the seed metal is lo- 
cated downstream of the arc discharge, it will not seriously alter the 
discharge operation. Therefore, the amount of metal seed that can be added 
will be nearly independent of the arc operation. Also, although the electron 
temperature in the arc exhaust is determined by the discharge conditions, it 
decreases with distance downstream from the active arc plume. As a result, 
by choosing the point at which the metal seed is added the electron tempera- 
ture that effects the gain can be controlled. An additional advantage of 
the MPD arc system is that it can be run at high powers (megawatts) so that 
large gains may be possible. 

An MPD arc system similar to that shown in figure 1 is being constructed. 
To estimate the possible gains that may be obtained a theoretical analysis 
was initiated. The following model was used to approximate the conditions 
that exist in figure 1. A beam of inert gas ions is assumed to be moving 
through a cloud of neutral ground state metal gas. Initially the ion beam 
and metal gas cloud have prescribed spatially uniform densities. Then 
treating this situation as an initial value problem, the densities of the 
important states are calculated as functions of time. The electron temper- 
ature is also calculated as a function of time. Conditions that exist in 
the initial penetration of the ion beam into the metal gas cloud should be 
approximated by this model. Results of this analysis will indicate the mag- 
nitude of the inversion, as well as the length of time the inversion will 
persist in the initial interaction of the ion beam and metal gas cloud. 
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III. ANALYSIS 

To obtain an exact solution for all the number densities of the ex- 
cited states of the metal ion would require solving the continuity, mo- 
mentum, and energy equations for each of the excited states of the metal. 
Also, similar equations for all states of the inert gas and electron gas 
must be included. Such a problem would be intractable. However, several 
simplifications can be made. First of all, only a limited number of levels 
of the metal ion will be important in determining an inversion. As a re- 
sult, most of the excited states can be neglected. Secondly, since elastic 
collision rates between like species are large, we can represent all states 
of the metal by a common temperature and velocity. Thirdly, the charac- 
teristic time for a change in density of an excited metal ion state is 
short ( 10 - 8 sec) compared to the momentum transfer collision time ( 10~6 sec) 
between the inert gas ions and the metal. As a result we assume that the 
velocity and temperature of the inert gas ions and the velocity and temper- 
ature of all the metal states are constants. Finally, heat conduction and 
viscous dissipation are neglected. With these approximations the system 
of equations to be solved includes a continuity equation for each consid- 
ered state of the metal, the electron continuity equation, the inert ion 
continuity equation, and an electron energy equation. 


til 

A. Continuity Equation for i Level of Metal Ion 

The continuity equation for the i^ level is given by the following 
expression (ref. 12, ch. 7-8). 


3n. 

i 

3t 


+ V • (n.v.) = P* + P 1 + P 1 , 
ii A e rad 


(3) 


where v. is the average velocity and n. is the number density of the 
ith species and t is time. The P terms on the right hand side of equa- 
tion (3) give the production rate of the i 1 -* 1 species as a result of charge 
exchange pf , collisional electron excitation Pg , and spontaneous 
emission and stimulated emission and absorption • In appendix A 

expressions for these terms are developed. In considering the charge ex- 
change production term we neglected the inverse deexcitation process (see 
eq. (la)). This is reasonable since the densities of the inert gas atom 
and excited metal ion are small compared to the inert gas ion and ground 
state metal atom. It was assumed that the radiation field used to calcu- 
late the stimulated emission and absorption effects is not dependent on 
the densities of the excited levels of the metal ion. Instead the radi- 
ation energy flux was assumed to have a Planck distribution in frequency 
and to be characterized by some radiation temperature, T r , and a param- 
eter, y, that accounts for an increase in energy flux at any frequency 
where lasing may occur (see appendix A). Since we have assumed v^ is 
nearly constant for time intervals of interest, equation (3) becomes the 
following. 
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3n. dn. . . . • 

I -> i i i i 

+ v. • Vn. = -7-^ = P^ + P + P . 
3t 1 l dt A e rad 


(4) 


B. Continuity Equation for Metal Ion Ground State 

For the ground state of the metal ion, an additional term must be 
added to equation (4) that accounts for ionization of neutral metal atoms 
by electron impact. Therefore, 


dn l _1 , _1 , .1 2,o 

-7— = P + P , + n n £ -n n,k , 

dt e rad o e eo e 1 el 


(5) 


Where subscript or superscript 1 denotes the ground state of the ion 
and subscript or superscript o denotes the ground state of the neutral 
metal atom. Note that i and are ionization and recombination 

rates. Similar to the case of electron excitation (eq. (A7)).we relate 
the recombination rate, k°^ , to the ionization rate, &g 0 , through the 
equilibrium relation (ref. 13, p. 687) 




2irm kTe 
e 


3/2 


qE 


I.P. 


kTe 


eo 


( 6 ) 


The quantity Ej.p. is the ionization potential of the metal atom in 
volts. Using equation (6) in equation (5) yields the following result. 


where 



(7) 


( 8 ) 


C. Continuity Equation for Metal Atom Ground State 

The continuity equation for the ground state of the metal atom in- 
cludes the effect of both charge exchange and electron collisionai ioni- 
zation. 



n 

e 1 



(9) 
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The first term in equation (9) is the sum of all the charge exchange 
rates (eq. (A9)) for levels of the metal ion that have a resonance with 
the inert gas ions. The second term is the net rate for electron col- 
lisional ionization. 


D. Continuity Equation for Ground State of Inert Gas 

The only loss mechanism for density considered for the inert gas 
ions was charge exchange. As a result the inert gas ion continuity equa- 
tion is the following. 

<10) 

r 


E. Charge Neutrality and Electron Energy Equation 

Finally, we assume that charge neutrality exists so that the elec- 
tron number density must balance the total ion density. Therefore, 


N 


n = 


n A + 


IK 

i=l 


(ID 


To complete the system of equations an electron energy equation is 
required. In Chapter 7 of reference 12 flow equations for general reac- 
tive systems are developed based on the Bhatnagar-Gross-Krook model for 
the Boltzmann equation. Based on these results the following electron 
energy equation was developed in appendix B from equation (56.2) of ref- 
erence 12. 



The term gg^ is the average of the square of the relative velocity be- 
tween all electrons and all level i . particles that take part in excit- 

ing level j. Similarly, gg^ is the average square velocity between all 
electrons and all ground state metal atoms that take part in ionization. 
The excitation rate of level j as a result of electron collisions with 
level i is denoted by . . ' 
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The system of equations is now complete. It consists of N - 1 
equations for the number density of considered excited levels of the . 
metal ion (eq. (4)), continuity equations for the metal ion ground state 
number density (eq. (7)), metal atom ground state number density 
(eq. (9)), inert gas ion number density (eq. (10)), and electron number 
density (eq. (11)), and finally of the electron temperature (eq. (B-5)). 
Before a solution can be obtained, however, expressions for the electron 
excitation rate, £^, electron ionizat ion rate, , charge exchange rate, 

£^ 0 , and the relative velocity terms g^. and g| must be known. In 


the next section expressions for these quantities will be developed. 


IV. CROSS SECTIONS AND REACTION RATES 


In appendix C the general expressions for the electron excitation 
and ionization rates (eqs. (C— 7 ) and (C-8)), are developed. The results 
are in terms of an integral over electron energy, <J> e , with the cross 
section for excitation, S3. , or -ionization, Sj, appearing in the inte- 
grand. To carry out the <j> e integrations in equations (C-7) and (C-8) 
expressions for the cross sections S^ . and Sj must be known. For 
ionization by electron collisions of flie metals to be considered (Mg, Ca, 
Sr, and Ba) experimental values of the cross section as functions of elec- 
tron energy are given in reference 14. From threshhold (cf> e = Ej.p ), the 
cross sections increase rapidly to some maximum and then begin to decrease 
slowly as <j> e increases. Since we will be considering electron tempera- 
tures below the value of <j) e for maximum cross section, a linear approxi- 
mation to the initial portion of the S-j- versus 4> e curve was used to 

calculate i?~ . 

eo 


; = a 1 (<f> - E_ D ) 

I eo T e I.P. 


*e± E I.P. 


S I = ° 


^e — E I.P. 


( 12 ) 


The quantity a is the slope of the initial portion of the Sj versus 
<|> e curve. Substituting equation (12) in (C-8) and integrating results 
in the following 



(13) 


In the case of electron excitation the Bohr-Thomson approximation 
for the cross section was utilized. Experimental values for these cross 
sections are not available. The Bohr-Thomson approximation used is of 
the following form (ref. 15). 
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(E. - E.) < 4 < (E , - E.) 
3 i ~ e - j+1 l' 


± ' E j + 1 - V 


(14) 


The parameter a Q = 6.513*10 in equation (14). Using equa- 

tion (14) in (C-7) yields the following result. 



The term Ei(x) is the exponential integral. 


Ei(x) 



(16) 


Equations (13) and (15) give the electron inelastic collision rates 
utilized in the solution. Since the Bohr-Thomson cross section for elec- 
tron excitation is usually larger than the actual cross section (ref. 15), 
we are overestimating the effect of electron collisions in the production 
of excited states of the metal ion. Similarly, since our approximation to 
the ionization cross section is larger than the experimental value we are 
overestimating the production of ground state metal ions. Overestimating 
the electron excitation and ionization rates underestimates the possibility 
of a population inversion by charge exchange. This is so because the elec- 
tron processes drive the metal ion le/el densities toward their equilibrium 
values thus eliminating the inversion. 


The calculation of the charge exchange rate, *Ao , is carried out in 
appendix D. In making the calculation a linear approximation for the 
cross section was utilized. 



Ao 


+ b I 

Ao 


( 17 ) 


are constants and <J> is the relative energy. 


Where af 
Ao 


and b^ 


9 


q<f> = t 


m.m 
A o 


2 m, + m 

A o 


(18) 


Values, for a\ n and b^ were obtained by fitting equation (17) to the 


_ Ao 0 ' — - 

data of references 16 and 17. Sifice the crata is for perfect resonance 

(A = B, AE = 0 in eq. (1)), we *&11 be overestimating the charge exchange 

rate. From equation (D- 6 ) of appendix B and using the approximation 


kT 


m. 


kT 


>> 


m 


we obtain the following expression for £ 


Ao‘ 


kT .m 

£ r = / A o r 

Ao V \q (m + m. ) a Ao 


3 (s 


c ,23 -2 

erf s| s + 4 " s + 


3 ) *!=-{'* I* 1 ) 


+ b 


Ao 


erf s 




-1 


/it 


(19) 


Where, 


-> 

v - v 
A o 


s - v _ | A o 
/2kTj /2kT^ 

\l^f V7T 


( 20 ) 


and erf s is the error function defined by equation (D- 8 ) . The quantity 
v will be referred to as the slip velocity since it indicates how fast 
the inert gas ion beam is "slipping" through the metal atom gas. 


The average relative velocity terms g|^ and g| Q are calculated 
using equations (B-3a) and (B-3b). Substituting equations (13) and (15) 
in equation (B-3a) and performing the integration yields the following 


m 


¥ - » . ■ f « (, j - v < ( Vi - v 


3k g ei 3 T e 9 l 


L 


kT, 


kT c 


E. - E. 

2 j+l 


=y 


(21) 


Where, 


8 , (x,y,z) ■ t . 2 . + -£ < 1 - n+y - 

e X - ze ^ - x[Ei(x) - Ei(y)] 


( 22 ) 


and Ei is the exponential integral defined by equation (16). Also, if 
equations (12) and (13) are substituted in equation (B-3b) the following 
is obtained. 
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m e 2 2 

3k 8 eo=3 T e 6 2 


qE 


I.P. 


kTe 


( 23 ) 


Where, 

2 

e 2 (x) — (24) 

Equations (21) and (23) were used in equation (B-5) to complete the elec- 
tron energy equation. 


V. MODEL FOR METAL ION 


There are many possible metal-inert gas systems where charge exchange 
may produce an inversion. However, the only system that is numerically 
evaluated in this study is calcium-xenon. A model of the calcium ion used 
in the analysis is shown in figure 2. Wavelengths for all the allowed 
transitions are given in nanometers and were obtained from reference 18. 
Only the 5s level has a near resonance with xenon ions [AE = -0.45 eV 
in eq. (1)]. As a result charge exchange will populate only the 5s 
level. This should be an advantage over other systems (refs. 1 to 8) 
where several Revels have near resonances and thus will compete for charge 
exchange population. 

As can be seen in figure 2, the two 4p, 5p, and 4d levels are 
nearly the same in energy. Using equations (A-3) and (A-6) we can write 
the net electron collisional excitation rate between two of these adja- 
cent levels as follows , 

q( VV\ 

kTe J ) > i (25) 

Referring to equation (15) we see that if E. -*E- then £*.-*“>. There- 
fore, in order that the term defined in equation "(25) remains finite the 
following condition must be satisfied. 


P J . . = n £ J . 
ei e ei 


in. - n . 

i 3 


8 i 


Therefore, 


g. 


q(E 1 -E ± ) 


n:, - 


n . 
3 


kTe 


-* 0 




q(E£-E.) 


'I kTe 
n. = n. — e 

1 J 8 i 


= n . 


£i 

j gj 


(26) 


(27) 
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which is just the Boltzmann equilibrium relation (eq. (A-6)). For the 
adjacent 4p, 5p, and 4d levels equation (27) was used to relate the 
densities. As a result the system of equations to be solved consists of 
N - L differential equations for calcium ion level densities, a differ- 
ential equation for xenon ion density, a differential equation for cal- 
cium atom density, a differential equation for electron temperature, and 
L algebraic equations for calcium ion level densities. The total number 
of calcium levels is N = 8, and L = 3 is the number of pairs of ad- 
jacent levels where the approximation given by equation (27) is used. 

The transition probabilities for the allowed transitions obtained 
from reference 19 are presented in table'll. Also appearing in table II 
are the ionization cross section coefficient, a^ 0 , and the charge ex- 
change cross section coefficients, a^ Q and b^ Q . The superscript 4 
refers to the 5s upper laser level. These coefficients were obtained 
by fitting the data of references 14 and 17. We can expect to find in- 
versions between the 5s level and the two 4p. levels. Inversions on 
these levels would produce lasing at 373 .7 nm and 370.6 nm. 


VI. RESULTS 

A numerical solution of the equations was carried out on an IBM 
7094 computer. Numerical approximations to the exponential integral, Ei, 
appearing in the excitation rate expression (eq. (15)) were obtained from 
reference 20. 

There are two important parameters that result from the calculations. 
The first is the gain defined as follows. 



(28) 


The subscript u denotes the upper laser level and £ denotes the lower 
laser level. For G > 0 an inversion exists and the gain coefficient, 
a, can be calculated using equation (28). For a Doppler broadened line 
the gain coefficient is given as follows (ref. 21), 


or , in MKS units , 


a 




g aVg 

u u 


a = 1.741x10 


i - _1 

h~ g A X G m 
1 T-. °u u 
15 


(29) 


The term m is the mass of an atom of the lasing gas, Tg is the temper- 
ature of the lasing gas, g^ is the degeneracy of the upper laser level, 
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is the transition probability for the laser transition, X is the 
wavelength of the laser transition, and Mg is the atomic weight of the 
lasing gas. The value of a for 10 ^ <. G <. 10^6 m -3 an( j x B = I.eV 
in the calcium-xenon system with X = 3737A° is 0.15 <. a<. 1.5 m - -*- 
and with X = 3706A°, 0.08_<a_<0.8. 


The second important parameter resulting from the analysis is the 
inversion time,, t. This has been arbitrarily defined as the length of 
time the gain, G, remains greater than lO^m - -^. The investigation con- 
sidered the effect of the initial conditions, (T e ) t=0 and (n 0 ) t _ 0 , as 
well as slip velocity, v, and charge exchange cross section coefficient, 
b^ Q , on the gain, G, and inversion time, t. For all cases investigated 
so far the xenon ion temperature has been arbitrarily set at T Xe - 0.1 
(T e ) t _ 0 . Also, the effect of the radiation field has been neglected, 

(T r -»■ 0). Finally, all densities except calcium atom density, n 0 , and 
xenon ion density (n x e)i> are zero at t = 0 . 


As mentioned before, inversions can occur between the 5s level and 
the two 4p levels. Since the densities of the two 4p lower laser levels 
are related by equation (27) the gain, G, and inversion time, t, for each 
of the laser lines will be equal (see eq. (28)). 


Figure 3 shows G as a function of time for three different ini- 
tial electron temperatures. The gain increases to a fairly "flat" maxi- 
mum and then decreases. For (T e ) t=0 < 3 eV a plateau occurs in the 
gain vs time curve. . However, for large electron temperatures this pla- 
teau disappears . We see that the inversion time is sensitive to initial 
electron temperature while the maximum gain is almost independent of 
(T e ) t _ 0 . These results are as expected. It is the electron temperature 
that drives the level densities toward their equilibrium values by elec- 
tron collisional excitation and thus destroys the inversion. Thus, the 
inversion time should depend strongly on the electron temperature. The 
maximum gain, however, is directly proportional to the maximum upper 
laser level density. Since this density depends mainly on the charge 
exchange mechanism, which is independent of T e , we do not expect the 
maximum gain to depend strongly on T e . 

The sensitivity of inversion time, r, to initial electron tempera- 
ture is shown in figure 4. We also see that t is much less sensitive 
to the initial calcium atom density, (n 0 )t=oj than to electron tempera- 
ture. An electron temperature, (T e ) t = 0 =0.7 eV, was the smallest value 
used in the calculations. Although x appears to be increasing as 
T e -*■ 0, we expect this trend to reverse for very small electron temper- 
atures. As T e •* 0, the deexcitation rate, k e x -> °°. Therefore, the 
upper laser level will be collisionally deexcited so fast that an inver- 
sion will not have time to form and t . -*■ 0. A complete explanation for 
the minimum that occurs near (T e ) t _ 0 = 1.5 eV has not been found as yet. 
This minimum disappears as (n 0 ) t=0 is decreased. The inversion time ! 
decreases rapidly from (T e ) t=:o =? 2 eV to (T e ) t=0 ^ 5 eV. Then the de^ 
crease is much slower and appears to approach an asymptote of about 10 "^ 
second. Figure 5 shows the dependency of x on the initial calcium 
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atom density, (n 0 ) t=0 . An increase of 10-*-9 to lO^l in (n 0 ) t ^ 0 decreases 
the inversion time by a factor of approximately 10. 


The maximum gain, Gu^, depends strongly on the initial calcium atom 
density, (n 0 ) t=0 . Figure 6 shows this dependence. As mentioned earlier, 
the maximum gain depends on .the charge exchange mechanism. From equation 
(A-9) we see that the charge exchange production term depends linearly on 
n o» n xe> and ^Ao‘ a resu l t we expect a strong dependence of G ma v 

on these quantities. Sinc e , the slip velocity, v, is greater than the 
xenon ion thermal speed, /2kTxe/ m Xe > the charge exchange rate, , can. 
be approximated by equation (D-9) . Therefore, is directly propor- 

tional to v. And since a^ Q is small, H^ Q is also directly propor- 
tional to 'bjL. Therefore, changing either (n 0 ) t=0 , (n xe ) t=0 , v, or 
b io win produce the same effect on the maximum gain; Figure 7 demon- 
strates this conclusion. A change of a factor of 10 in either (n 0 ) t=0 
or v produces nearly the same change in G,,^. . A plot of vs 

h^o or ^max vs ^ n xe^t=o wou ld show results similar to those of fig- 
ure 6 or 7. 


Although not shown in this report, the analysis predicted the fol- 
lowing results for electron temperature. For (T e ) t _ 0 < 5eV the elec- 
tron temperature remains nearly constant with time. At (T e ) t=o = 5eV 
the temperature decreases by less than 5 percent for the time interval 
0 <. t <. t. For larger initial electron temperatures the decrease becomes 
significant. When (T e ) t=0 = 20 eV the electron temperature decreases by 
25 percent for the time interval 0 <_ t <. x. 


VI. CONCLUDING REMARKS 


Results obtained for the calcium-xenon charge exchange system indi- 
cate maximum gains of lO-^-lO^ m - ^ with inversion times of 2x10"® - 
8x10~7 second. In obtaining these results the effect of the radiation 
field. (stimulated emission) has been ignored. Including the radiation 
will certainly reduce the gain. 

The inversion time, x, is determined mainly by the electron colli- 
sional processes and thus by the electron temperature. For electron 
temperatures greater than.2eV, x decreases rapidly. Maximum gain, G max , 
however, is nearly independent of electron . temperature and depends mainly 
on the charge exchange mechanism. . As a result, increasing either the in- 
itial xenon ion density, initial calcium atom density, slip, velocity, or 
charge exchange cross section will yield the same increase in the maxi- 
mum gain. 

As pointed out previously, we are overestimating the charge exchange 
rate by. using the resonant cross section to approximate the near resonant 
case. Therefore, the results for maximum .gain are optimistic. Also, 
however, the Bohr-r Thomson cross section used overestimates the electron 
excitation processes. As a result . the inversion times calculated in the 
analysis are underestimates of actual inversion times. 
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Consider how the results apply to the MPD arc system described in 
figure 1. As mentioned in the Introduction, the model used should ap- 
proximate only the initial interaction between the xenon ion beam and 
the neutral calcium gas cloud. An estimate of the length of the initial 
interaction for G > 10^ m“3 is given by the product of the inversion 
time, T, and the slip velocity, v. Therefore, from the results we ex- 
pect the initial interaction to have lengths the order of 0.05 - 1 cm. 
The inversion will not necessarily disappear after this initial inter- 
action in the actual system, however. Since the xenon ion beam will 
continue to penetrate into the calcium cloud, we would expect the inver- 
sion to be maintained until the calcium cloud is completely penetrated 
by the xenon ion beam. The gain, however, will be less than the gains 
calculated for the initial interaction region. 
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APPENDIX A - DEVELOPMENT OF PRODUCTION TERMS FOR iTH 
LEVEL CONTINUITY EQUATION 

1. Collisional Electron Excitation 


Population of the 
and radiative effects, 
for exciting the ith 


ith level occurs as a result of both collisional 
We consider two possible collisional mechanisms 
level; electron collisions 


e + B. 


ej 

k r 

ei 


e + B. 
1 


(A-l) 


and if the level is near resonant then charge exchange described by equa- 
tion (la) is also important. In equation (A-l) £|. is the electron col- 
lisional excitation rate of the ith level from tne jth level and 
k3 ei is the corresponding deexcitation rate from i to j , (i > j). 

The contribution to the density of the ith level from the reaction 
described by equation (A-l) is the following (ref. 12, p. 264). 


n nj 1 , - n n,k , i > j (A-2) 

e j ej e i ei 


where n denotes density and the subscript e denotes electrons. Term 
(A-2) gives the contribution to the ith level from lower levels (i > j). 
Similarly, the ith level contributes to the density of all levels above 
the ith. 

- (n ri . £^ . - n n . k 3- j > i (A-3) 

\ e l ex e l ejy J 

Term (A-3) gives the loss rate of the ith level due to excitation of a 
level j above the. ith level, (j > i) . Therefore, the net production 
of the ith level due to electron excitation is the following. 



i-1 




N 


n e n i k ei) ’ D 


n n . £ . 
e l ei 


j=i+l 


n n.k^.j . (A-4) 

e J e 3/ 


Where N is the total number of excited levels considered. (The ground 
state of the metal ion is denoted by i = 1.) 

If equilibrium conditions exist then the production rate terms must 
vanish. For example, equation (A-2!) set equal to zero yields the 
following 
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(A-5) 


At equilibrium the densities obey the Boltzmann relation (ref. 13, p. 63) 




q(E i -E.) 

kT 

e 


(A-6) 


Where g denotes degeneracy, E is the energy of a given level in eV, q 
is the magnitude of the electron charge (1.6x10^-^ coulombs), and k is 
the Boltzmann constant (1.38x10“^^ joule/°K) . If we substitute equation 
(At 6) in equation (A-5) the following is obtained. 



eq. 




q(E.— E . ) 
.3 i 

kT 

e 


(A-7) 


Using the usual assumption that equation. (At 7) holds away from equilib- 
rium we can replace k£^ and k*j in equation (A-4) by equation (A-7). 
Therefore we get the following. 

qE i qE, 

i-1 


P 1 = n 
e e 


kT i-1 


?n/ \ P^.g. 

Lj j ej g, i L/ ej J 

j=l j=l 


kT 


q E i 

n j kT e V\ j 

- n i b l ei + g ± e b t *ei € 


qE 

kT 


j=i+l 


j=i+l 


(A-8) 


Equation (A-8) gives the production rate for the ith level as a result 
of collisional electron excitation. 


2. Charge Exchange Excitation 

If the ith level is in resonance or near resonance with the inert 
gas ion, then charge exchange will also contribute to the production of 
the ith level. 
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n.n Z , 6 . 
A o Ao rx . 


(A-9) 


Where is the density of the inert gas ions and n Q is the density 

of the ground state of the metal atom. The subscript r denotes a level 

in resonance or near resonance with the inert gas and 6 r ^ = 0 if r 4 i 

or 6 . =® 1 if r = i. 
rx 

3. Spontaneous and Stimulated Emission 

Radiation effects contribute to the density of the ith level by 
way of spontaneous emission and stimulated emission and absorption. The 
net contribution from radiation is given as follows (ref . 12, p. 320) 


i-1 

’rad ' C h B ) ^ - “i A i] 

j-1. 


IN . 

+ n i B i^ + n j B j<t> ( v |) + n j A ^J (A-10) 


j=i+l 


The term A | is the spontaneous emission coefficient (or transition 
probability) with the units of sec”! for spontaneous emission from the 
ith level to the jth level. The stimulated emission or absorption co- 
efficient for stimulated emission or absorption from the ith level to 
the jth level is denoted by B^. It has units m^/joule sec. The ra- 
diation field energy flux per frequency interval (joule/m^) is given by 
<p (vj) where 


hv^ = q(Ej 


E.) = q E. - E. 
x n i j 


(A-ll) 


The first sum in equation (A^IO) gives the radiation contribution to the 
ith level from all levels below i. The first term in the sum being 
stimulated absorption- where a lower, level, is. converted to the ith level, 
the second term is stimulated, emission where the ith level drops to a 
lower level, and the last term is spontaneous "emission where the ith 
level drops to a lower level. The last sum in equation (A-10) gives the 
radiation contribution to the .. ith level, from all levels above i. Each 
term has a meaning similar to those in. the first sum except that i now 
replaces j as the lower level. 

For an accurate solution to the problem the radiation field energy 
flux, 4> (v), should be coupled to the densities of all species through a 
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radiative transfer equation.. However, to simplify the problem the fol- 
lowing form for <j> (v) was assumed. 


<j> (vj) 



-1 


(1 + 


yd . 6 . . ) 
uj £i' 


(A-12) 


In equation (A-12) c is the speed of light and <Sfj = 0 for i ^ j 
and 6^j = 1 for i = j. The subscripts u and £ denote levels where 
a population inversion exists . Thus, u represents the upper level and 
£ the lower level of the inversion- In choosing equation (A-12) for 
<j> (v) we are assuming the radiation field has a Planck distribution ref. 

13, p. 721) in frequency and can be characterized by some radiation tem- 
perature, T r (ref. 12, p. 320). In addition, the Planck distribution has 
been multiplied by the factor (1 + yd u jd£i) where y is a constant. This 
factor is added to account for an increase in radiation energy density 
that will result if lasing occurs on. the transition, u £. Both y and 
T r are treated as constants in the analysis. 


To simplify equation (A-10) we use equation (A-12) for <f> (v) and the 
Einstein relations (ref. 13, p. 722), 



g i 1 




2hv“! 

i 


— A 1 
3 


(j > i) 


(A-13) 


where g^ is the degeneracy of the ith level. Using equations (A-ll) , 
(A-12), and (A-?13) in. equation (At- 10) yields the following result. 


N 


N 


rad 




j=i+l 


n.A. (1 + yd .6 ) 

.1 .1 uy £i 

exp[q(E i - E ± )/kT r ] - 1 


i-1 


j=i+l 


C n.A. (1 + yd .d ) i r — v 

]. l . m £] _ J \ ' 

g (exptq(E i - Ej ) /kT r ] - 1} n i| j_j 

j=l J 3=1 


i-1 


i-1 


N 


+ 


C Aj(l + yd .5 ) g .At (1 + yd . d ) i 

_x_ ui £j _L ) ' 3 3 uj £i l 

exp[q(E i - Ej)/kT r ] - 1 g . l_J exp[q(E i - Ei)/kT r ] - 1[ 

j=l j=i+l 1 


(A-14) 

Equation (A-14) gives the production rate of the ith level by spontan- 
eous emission and stimulated emission and absorption. 
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APPENDIX B - ELECTRON ENERGY EQUATION 


From equation (56.2) of reference 12 we obtain the following electron 
energy equation. 


dT l 

n k — — = — n k 
e dt 2 e 


3T 

P — 

— - + v • VT 
3t e e 




K 

es 


3k (T - T ) 
s e 


m 


N N-l 


D E 

j=i+l i=l 


q(E.-E.) 

liJ. 1 


8 i 

n. - — n .e 

1 8 j J 


kT 


m 


T 4 - I kT e 


n £ 
e eo 


n - n n.,Y , 
o e l /v, el 


e 2 3 

T^eo-I kT e 


(B-l) 


Where, the first term on the right-hand side of equation (B-l) is the 
elastic electron collisional energy transfer rate. The second term is 
the energy transfer rate due to electron excitation collisions of the 
metal ion, (eq. (A-l)). And the third term is the energy transfer rate 
resulting from electron collisional ionization of neutral metal atoms. 
Other terms appearing in equation (B-l) are the following. 


8kT 


K = — m n n 
es 3 e e s \/ ^nn 0 es 


(B-2a) 


CT es 2 IkT 



q^ e 

0 kT 

m 2 e,, 

a d> e d <f> • 

es e e 
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2 kT 
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i e e 


(B-3a) 


(B-3b) 
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Definitions of these terms are obtained from reference 12. Three cross 
sections appear in equations (B-2b) and (B-3b) . They are the elastic 
momentum transfer cross section, o® g , for collisions between electrons 
and species s, the inelastic excitation cross section, . , for produc- 
tion of level j due to collisions between electrons and e level i, and 
the electron ionization cross section, Sj, for the metal atom. All of 

these quantities are functions of electron energy, <fu . The term, 
i i 

g|^, is the average of the square of the relative velocity between the 
electrons and level i particles that take part in excitation of level j 

Also, gg Q is a similar average square velocity for electrons and neutral 
metal atoms that take part in ionization. At equilibrium the following 
condition holds for a reaction like that described by equation (A-l) . 


“is 2 
2 g ei 


m 


LJ 
e 2. 

2 8 


+ (E. - E.) 
ej J i 


(B-4a) 


Also , at equilibrium 



U 

2 


m '-z 

"T 8 el +E !.P. 


(B-4b) 


The quantity (E-s - E^) is the energy required to produce level j by 

excitation of level i and is the average of the square of the 

relative velocity between, the electrons and level j that result from 
the excitation reactions. The ionization potential of the metal atom 

is ^I.P. and b|i i s the average square relative velocity between 
electrons and ground state metal ions that result from ionization. In 
obtaining equations (B-l) through (B-3) the following approximations 
were made. 


1. m + in 
e t 


m. 


kT kT 

— ® t 


m 


in 


t 


2. kT >> m^v 
e t t 


(v is average velocity of species t) 


3. Equations (BT-4a) and (B-4b) 


4. Heat conduction negligible 

5. Velocity gradients negligible 


Referring to equation (B-2a) we see that the elastic transfer rates 
(first term on right-hand side) will be. proportional to the mass ratio, 
m e /m g . Since the mass ratio is small the elastic energy transfer rate 
was neglected. As a result equation (B-l) becomes the following 
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APPENDIX C - ELECTRON EXCITATION AND. IONIZATION RATES 


A general expression for the production rate of species r as a re- 
sult of collisions between species s and t based on the Bhatnagar- 
Gross-Krook model for the Boltzmann equation is given in reference 12 
(p. 267-268), 

2 
v 
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st 3/2 3 
it a- 


oo 
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r j 

n 2 tt 

' g 3 Sg t (g)exp 

2 

g 2gv cos 0 

2 2 
a a 

OQ 

II 

O 

C 

9=0 U 

o 

ii 

■©- 

J 


sin 9 d0 d <f>dg 


(C-l) 


Terms appearing in equation (C-l) are the following 

- relative average velocity 


v = v - v 
1 t s' 


(C-2) 


g = l c t - c s l 


- relative velocity between colliding 

pair of particles (C-3) 


„ 2kT 2kT 

2 S , t 

a = k 


m 


in 


(C-4) 


The cross section for the reaction, S s t, is a function of the relative 
velocity, g, only. As a result the angular integrations in equation 
(C-l) can be performed. 
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,vtT 



2gv _ 2gv\ 
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a< 2 „r 

6 S st d 8 


(C-5) 


Consider equation (C-5) for electron excitation or ionization col- 
lisions. In this case, 


kT kT 

— ® t 


m 




and 
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v 

a 



<< 


1 


m 


As a result the exponential terms in equation (C-5) can be approximated 
by 


2gv 


e 


2 


a 




(C-6) 


Using equation (C-6) in (C-5) and also changing variables from g to 

12 

electron energy q<j> £ = y m^g yields the following results for excita- 
tion and ionization. 



(C-7) 


(C-8) 
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APPENDIX. D - CHARGE EXCHANGE RATE 

We begin the calculation of the charge exchange rate by rewriting 
the general rate equation (eq. (C-5)) in the following manner. 


Ao 


olv/tt 



_ (v+g)' 


2 

a r 

S Ao d S 


1 v o 


(D-l) 


Now make the following variable changes.- In the first integral let 
y = v - g and in the second integral let y = v + g. Also, use the 
condition that S][ 0 (g) = Sj[ 0 (-g), since it is a function of relative 


energy, T 


- m.m 
1 .. A.o 


2 m. + m 
A o 






(D-2) 


To proceed further an. expression for the cross section, S^ Q must 
be known. As pointed out in the Cross Sections and Reaction Rates sec- 
tion, a linear approximation was used for S£ 0 . From equations (17) 
and (18) 



r 1 2 

a Ao 


+ b 


r 

Ao 


where the reduced mass is given by 
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U = 


m.m 
A o 

m. + m 
A o 


(D-4) 


In obtaining equation (D-2) we have used a change of variables such that 
g2 = y2 _ 2y v = v^. As a result, the cross section to be used in equa- 
tion (D-2) has the following form. 


„r r u , 2 „ , 2. , ,r 

S Ao ' a A0 <y ‘ 2vy + V > + b * 


Ao 


(D-5) 


Substituting equation . (D-5) in (D-2) and integrating yields the 
following result. 
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and 


erf S = 


2 

-x 

e dx 


(D-8) 


It is interesting to look at two limiting cases for equation (D-6). When 
the difference in average velocities, v, is large compared to the average 
thermal speed, a, then S -> 00 and equation (D-b) becomes the following. 
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(D-9) 


When a » v, equation (D-6) becomes the following. 



For the case S the charge exchange rate is directly proportional to 

the difference in average velocities. When S .-*• 0, the charge exchange 
rate is directly proportional to the average thermal speed. 
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TABLE I - POSSIBLE CHARGE EXCHANGE LASER SYSTEMS 


Upper Laser Level . 
Pumping Reaction 

Laser 

Transition 

Laser 

Wavelength , 

X, 

nm 

Energy 
defect , 
AE , 
eV 

Efficiency, 

hc/XI.P. 

Xe + + Ca Xe + (Ca + ) u 

5s 2s i /2 - 4p 2 P 1/2 

.370.6 

-0.451 

0.27 

5s 2 S;l /2 - 4p 2 P 3/2 

373.7 

-0.451 

0.27 


3p2p i/2 * 3 S 2s i/2 

280.3 

0.0615 

0.36 

Xe + Mg -*■ Xe + (Mg ) 

3p2p 3/2 + 3s2s i/2 




u 

279.5 

0.0501 

0.36 


5d 2D 5/2 -► 5 P 2p 3 /2 

346.4 

-0.182 

0.29 

Xe + + Sr -> Xe + (Sr + ) 

u 

5d - 2 ' D 3/2 ** 5p2P 3/2 

347.5 

-0.17 

0.29 


5d2D 3/2 + 5p2p i/2 

338.1 

-0.17 

0.30 


5 £2 F 5/2 - 6d 2 D 5/2 

889.6 

-0.197 

0.11 

Xe + + Ba -> Xe + (Ba + ) 

u 

5£2f 3/2 - 6d 2 D 3/2 

873.4 

• 

-0.197 

0.11 


5£ 2 F 7/2 - 6d 2 D 5/2 

870.9 

-0.227 

0.11 

Ar + +Mg -*• Ar + (Mg + ) 

4s 2 Si /2 * 3p 2 P 1/2 

292.9 

-0.542 

0.27 

4s2s 1/2 3p2P 3/2 




u 

293.6 

-0.542 

0.27 

+ 4- 

4£ 2p 7 / 2 * 4d2D 5/2 

892.7 

-0.002 

0.09 

N + Ca -* N + (Ca ) 

u . 

4f2p 5/2 +4d2D 3/2 

891.2 

-0.002 

0.09 
































TABLE II - PARAMETERS FOR CALCIUM MODEL ION 



Charge exchange cross section coefficients: a. = - 2.98x10 
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A + + B -* A + (B + )j 

Figure 1. - MPD arc charge exchange laser system. 
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INVERSION TIME, x, SEC 



Figure 4. - Inversion time vs initial electron 
temperature for Ca-Xe. 
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Figure 5. - Inversion time vs initial calcium num- 
ber density for Ca-Xe. 



Figure 6. - Maximum gain vs initial calcium 
number density for Ca-Xe. 



SLIP VELOCITY, v, m/SEC 


Figure 7. - Maximum gain vs slip velocity for Ca-Xe. 



